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Abstract 

A mixture of spin- 1/2 fermionic atoms and molecules of paired fermionic atoms 
is studied in an optical lattice. The molecules are formed by an attractive nearest- 
neighbor interaction. A functional integral is constructed for this many-body system 
and analyzed in terms of a mean-field approximation and Gaussian fluctuations. This 
provides a phase diagram with the two merging Mott insulators and an intermediate 
superfluid. The Gaussian fluctuations give rise to an induced repulsive dimer-dimer 
interaction mediated by the unpaired fermions. The effect of an unbalanced distribution 
of spin-up and spin-down fermions is also discussed. 

1 Introduction 

A wide new field for investigating complex many-body systems has been opened by the 
idea that clouds of atoms can be cooled to very low temperatures by sophisticated cooling 
techniques pQ. Such an atomic cloud can be brought into a periodic potential which is created 
by counterpropagating laser fields [2]. This potential mimics the lattice of core atoms of a 
solid-state system and is called optical lattice due to its origin. It allows the simulation of 
conventional condensed-matter systems as well as the creation of new many-body systems. 
New quantum states can emerge due to the interplay of tunneling and interaction between 
the atoms. 

In the following, a cloud of spin- 1/2 fermionic atoms in an optical lattice is considered, 
where an attractive interaction between atoms in nearest-neighbor lattice wells is assumed. 
Our aim is to study different quantum phases that can appear due to the formation and 
dissociation of molecules, and condensation of the (bosonic) molecules. 



2 Model 

Our atomic cloud in an optical lattice is described by a grand-canonical system of spin-1/2 
fermions at temperature 1/(3 and chemical potentials /ii and fi2, referring to the two possible 
projections of the spin. Its partition function is defined as a functional integral with respect 
to Grassmann fields as [31 H]: 

Z = fe s ^D[M], (1) 
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where S is the action 



dr. 



(2) 
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t is the tunneling rate of single fermions, whereas J is the rate for tunneling of a pair of 
fermions, located at nearest-neighbor sites in the optical lattice. The J-term represents 
an attractive interaction. In contrast to a local interaction, it provides a dynamics for 
the molecules. Depending on the ratio t/J, there is a competition between the individual 
fermion dynamics (dominating for t/J^> 1) and the dynamics of molecules (dominating for 
£/J<l). 



3 Mean-field approximation and Gaussian fluctuations 

We decouple the J-term in the functional integral by two complex fields <fi and x [HI Hj- 
cf) is related to the order parameter for the formation of molecules, and \ is required for 
stabilizing the complex integral. A subsequent integration over the Grassmann fields leads 
to 

SeS = J |E ^rV'r'^r' + XrXr - In det G 1 > dr , (3) 
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and a nearest-neighbor matrix w whose elements are l/2d on the d-dimensional lattice. 
Moreover, we have v = J(w + 21). The saddle-point (SP) approximation 5S e s = for 
uniform fields and \i\ = /12 leads to the BCS-type mean-field result: 

X=-^, k = G, (4) 
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with 

G — — - / / ; .„ , -, t~, —dx , (5) 

P J ~ 1 l^l / 9 _ \ lUJ n + Hi + tx){lUJ n -fl2~ tx) 

and the density of states p of free particles in the optical lattice. The mean-field calculations 
gives three phases: an empty phase, a Mott insulator and a Bose-Einstein condensate (BEC) 
of the molecules whose condensate density is 

The latter is plotted in Fig. (TJ and the phase diagram is shown in Fig. [2j 

Excitations out of the molecular BEC can be described by Gaussian fluctuations around 
the SP solution by complex fields and x- Using A = i(p + x, A = i<fi + x, the corresponding 
action is 

<Seff — S'eff + fiSeft (7) 
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The above result reads in terms of Fourier coordinates 

SS ef£ = ^(6<i> q>ul ,G^(q,iu)6<f> q}U ), 



(9) 



q,u> 



where G^ 1 is a 4 by 4 propagator, 5<p is a four-component spinor. More details of the 
calculation can be found in Ref. pE]. The excitation spectrum e q = iu(q) is the solution of 



det G;£(q,iu) = 0. 



(10) 



This gives the Bogoliubov spectrum in the molecular BEC and a gapped spectrum outside 
the BEC (see Fig. EJ. 



4 Discussion and Conclusions 

Our mean-field approach also allows us to consider an unbalanced molecular condensate 
[5. 6 J by imposing different chemical potentials for the two spin projections of the fermions. 
A*i = A* + h, [1% — n — h. Although we cannot address questions about nonlocal properties 
(like phase separation [3, El [9]) directly within our mean- field approach, the effect of two 
chemical potentials provides interesting effects even in a uniform system. In particular, the 
existence of a first-order phase transition, usually leading to phase separation, can be studied 
with the mean-field action of the unbalanced system 



S cS ~ — --J_ iP (x)ln 



where 



dx, (11) 



E±(x) = -h± 710179 + (/i + tx) 2 . (12) 

This gives a first order-phase transition due to two separated minima in the mean-field action 
for small \x or for larger fi, depending on h and t (for fixed J) (cf. Fig. Hj). 

For small single- fermion tunneling rate t, there is a spin-polarized phase simply because 
one spin projection has a negative chemical potential. This can happen for small \i (cf. Fig. 
[5]). If t is larger, there is a shift of the chemical potential in Eq. ffl2l by the single-fermion 
tunneling rate. This prevents the appearence of a spin-polarized state for small [i but it 
leads to a sudden disappearence of the molecular condensate at large as shown in Fig. 
This is accompanied by a first-order phase transition. At the point h = \i there might be 
a coexistence of molecules and spin polarized fermions [TU1 [TT] . 

In conclusion, our mean-field approach to the model of Eq. (|2J) reduces to a BSC-type 
theory of molecules for spin-1/2 fermions with attractive interaction. It also provides a 
Mott insulator and a spin-polarized phase. The Gaussian fluctuations describe Bogoliubov- 
type excitations of a molecular condensate and the gapped excitation spectrum of a Mott 
insulator. 
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Figure 1: Density of a molecular condensate for different single-fermion tunneling rates t. // 
and i are in units of J. 




Figure 2: Phase diagram for different values of t. 
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Figure 3: Quasiparticle excitations for t = 0.5 and J = 1. 
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Figure 4: Mean-field action of an unbalanced system, indicating a first-order transition from 
the molecular condensate with increasing h = (fii — For t = 0.2, /x = and zero 

temperature the transition takes place around h = 0.25. 
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Figure 5: Density of a molecular condensate as a function of (/ii + //2)/2 for h = 0.26, J = 
1,T = 0. 
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